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Abstract. Previous work on the neutrino spectra from high 
energy sources has not included the tau neutrinos directly 
produced by the decays in the source. Here we consider the 
tau neutrino component and discuss how its inclusion mod- 
ifies the expected neutrino spectra. We discuss implications 
for interpreting any observed tau neutrino component in TeV 



- UHE events as evidence of v.. 



oscillations. 



1 Introduction 

The observation of high energy neutrinos from beyond Earth 
will open a new window on astrophysics and cosmology. 
Presently much effort is being expended on designing and 
deploying detectors and methodologies capable of resolving 
the flavour of the incoming neutrino flux. One motivation for 
doing this is that the absolute sensitivity, or relative sensitiv- 
ity compared with background, may be significantly greater 
for one neutrino species than the other species (Dutta et al., 
2000). Another primary motivation for identifying the neu- 
trino flavour is to test the hypothesis that neutrinos undergo 
flavour oscillation as they propagate from source to Earth. 
Indeed a number of atmospheric, accelerator and solar neu- 
trino experiments hint at the existence of neutrino oscilla- 
tions, although the experiments are not yet sufficiently con- 
sistent to delineate a unique solution (Akhmedov, 2000). The 
most favoured solution, fitting the latest SuperKamiokande 
data, involves — > v T oscillation with a mass difference of 
AmJ ~ 10 _3 eV 2 . Because the non-oscillation production 
of v T is assumed to be negligible relative to v e and produc- 
tion, and the cascade v r produced by collisions of primary v e 
and Vfj, flux with the relic cosmic neutrino background is or- 
ders of magnitude less than the primary v e and flux, any 
detected v r component from a source outside the Solar Sys- 
tem is currently expected to be indicative of oscillation. 

In this paper, we briefly report some results of our inves- 
tigation into high energy v T production, presented elsewhere 



(MacGibbon et al., 2001). We illustrate that the v T compo- 
nent at the high energy end in the hadronic decays is signifi- 
cantly higher than previously assumed and that this can have 
observational consequences. 



2 Neutrino Production 

A number of sources for TeV-UHE neutrinos have been pos- 
tulated: for example active galactic nuclei (AGN) (Athar et 
al., 2000; Husain, 2000); topological defects such as super- 
conducting, ordinary or VHS cosmic strings (Bhattacharjee 
and Sigl, 2000; Hill et al., 1987; Wichoski et al., 1998); su- 
permassive gauge and scalar particle (X particle) decay or 
annihilation (Albuquerque et al., 2000); and Hawking evapo- 
ration of primordial black holes (MacGibbon and Carr, 1 99 1 ; 
MacGibbon and Webber, 1990; Halzen et al., 1995). In all of 
these scenarios, the final state particle distributions are ex- 
pected to be dominated by hadronic decays at the source. 
In these decays over 90% of the final state emission will 
be 7T°, 7r + , and 7r _ , with the remainder mainly nucleons 
which decay into protons and antiprotons. On astrophysical 
timescales, the tt^ decay into v e and and e ± . It is known 
that the final cluster states (pions and nucleons) in QCD jets 
at accelerator energies are well described by the fragmenta- 
tion function 



™ = 15 ^ _ 
dx 16 



(1) 
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where x = E/mj and mj is the total energy of the decaying 
jet (Hill et al., 1987). This distribution continues down to 
E ~ 1 GeV. When convolved with the decay, Eq.(l) 
leads to, similarly, a dominant E~ 3 / 2 term in the v e and 
spectra at £ < 0.1 and dN VejV JdE — > as x — ► 1. (For a 
full derivation of the v e and spectra using fragmentation 
function (1) see (Wichoski et al, 1998).) 

The production channels for the v T , however, are substan- 
tially different. The tau neutrinos are only produced in signif- 
icant numbers once the tau lepton and heavy quark masses, 
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Fig. 1. The neutrino spectra, dN/dE vs x — E/mj, generated by 
the decay of raj = 10 TeV qq jets. The solid line represents v^, 
the dashed line represents v e , and the dotted line represents v T . 



Fig. 2. The E 2 dN/ dE vs x — E/mj neutrino spectra for the same 
jet decay as in Fig.(l). 



m T = 1.78 GeV, m b ~ 4 GeV, m t ~ 174 GeV, are sur- 
passed. The greatest contribution comes from the t quark 
decay. Because v T production is suppressed compared with 
other species, it has been assumed in flux calculations that 
the v T spectrum is orders of magnitude less than the v e and 
spectra at all x. This is not so. While indeed the total 
number of v T produced per jet is less than 10~ 3 of the total 
number of i/ e and v^, the high x tau neutrinos are predom- 
inantly produced by the initial decays of the heavier quarks 
with shorter lifetimes and the v e and are produced by the 
final state cluster decays of the much lighter pions. This leads 
to significantly greater relative contribution from v T at high 
x than previously assumed. The fragmentation distribution 
(1) is no longer relevant for the tau neutrino. 

In Fig.(l) we show the v T spectrum, together with the v & 
and spectra, generated by the decay of 10 TeV qq jets. 
To simulate these spectra we used the QCD event generator 
HERWIG (Corcella et al., 2001) and the process e+e" -> 
QiQi (<?)> « = all q flavours. Consistent spectra are obtained 
with PYTHIA/JETSET (see MacGibbon et al. (2001)). Note 
that in the region 0.1 < x < 1, the tau neutrinos make 
up more than one tenth of the total neutrino contribution. 
We also note that below x < 0.1, dN v JdE falls off with 
roughly an E~ x / 2 slope, and not the E^l 2 slope of the 
v e and spectra. Because the 1 GeV - UHE cosmic ray 
flux backgrounds are expected to fall off as dN/dE oc E~ y 
where 2 < y < 3, we also plot E 2 dN/dE in Fig.(2). 

From our 300 GeV - 75 TeV simulations, we find that the 
v T spectrum generated by qq jet decay can be parametrized 
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(2) 



per jet. Similarly the number (multiplicity) of neutrinos pro- 
duced per qq jet and the average neutrino energy scale as 
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Note that the average v T energy is significantly higher that 
the average v e and energies, consistent with our remarks 
above. We derive similar dN Vr / dE VT , N Ut and E Vt param- 
etrizations for the v T generated by initial tt decay, bb decay, 
t + t _ decay, W + W~ decay and decays which include ex- 
tensions to the Standard Model (e.g. SUSY and Higgs sec- 
tors) in (MacGibbon et al., 2001). 

In analogy with the use of the fragmentation function (1), 
we now apply (2), and its high energy extrapolation, to in- 
vestigate the neutrino fluxes produced in astrophysical and 
cosmological scenarios. 



3 Neutrino Fluxes from Astrophysical and Cosmologi- 
cal Sources 

To calculate the expected v T flux from a given astrophysical 
or cosmological source, the above v r fragmentation function, 
Eq.(2), must be convolved with the initial distribution of the 
decaying particles, the Galactic or cosmological distribution 
of the sources and the redshift evolution and interactions of 
the emitted neutrinos as they propagate from source to Earth. 

Here we briefly mention two scenarios. The new spectra 
for the tau neutrino fluxes expected from a number of other 
astrophysical and cosmological sources, including AGN and 
primordial black holes, are presented in our paper (MacGib- 
bon et al.,2001). 
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Fig. 3. The neutrino spectra from VHS cosmic strings with a mass 
per unit length of G/i = 10~ 8 . The dash-dotted curve is the v T 
component produced by cascades off the relic neutrino background 
(Yoshida et al., 1997). 
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Fig. 4. The neutrino spectra from VHS cosmic strings with G/Lt = 
1CP 10 . Points with arrows represent upper limits for the flux (taken 
from (Bhattacharjee and Sigl, 2000)). 



Our first example is the VHS cosmic string scenario (Wi- 
choski et al., 1998). However, our remarks apply generically 
to p = 1 cosmic string models. In the VHS scenario, 1 GeV - 
UHE particle fluxes are generated by the decay of supermas- 
sive scalar and gauge particles emitted by the long strings 
over the age of the Universe. The neutrino spectra expected 
at Earth from VHS strings with a mass per unit length of 
G/j, = 1(T 8 are shown in Fig.(3) (Wichoski et al., 1998). 
The fraction of the v T component which is solely produced 
in the collisions of the primary v e and with the relic cos- 
mic 1.9° K neutrino background is represented by the dash- 
dotted line. (This is the only v T component presented in pre- 
vious cosmic string papers.) The dotted line represents the v r 
directly produced in the hadronic decays of the string emis- 
sion. As can be seen, our results give a significant increase 
in the expected v T signal at the highest energies. Note too 
that because the cosmic ray background falls off as E~ v , the 
highest energy region of the spectra is most relevant to detec- 
tion. In Fig. (4), we show the neutrino spectra from a viable 
VHS model with G/i = 10~ 10 and the comparisons with 
present and proposed detector sensitivities. 

As the second example, we briefly comment on the im- 
plications of our results for the simpzilla scenario. The an- 
nihilation of strongly interacting superheavy dark matter X 
particles captured by our Sun (simpzillas) has recently been 
explored as a source of observable v T (Albuquerque et al., 
2000). In this scenario the X particle decays into 2 quarks 
or 2 gluons which then decay into jets. Because of the high 
solar density, it is argued that only the t — > W — > v chain 
in the X decay would produce a significant v T flux capable 
of escaping the Sun. Assuming the fragmentation function 
(1) applies to the t quarks in this decay chain, the authors de- 
rive a distribution for the v T produced per X particle decay 
which has a leading E~^l 2 slope above m t . However, our 
full W decay simulations show that the v T spectrum falls 
off with a much weaker slope, inconsistent with an E~ 3 / 2 



slope from the partial decay. As a resolution to this discrep- 
ancy, we believe that the application of Eq.(l) to the initial 
t distribution from X annihilation in the Sun is inappropri- 
ate because the short-lived t quarks are generated in the first 
step directly by the X annihilation, whereas Eq.(l) is a fit 
to the final state cluster distributions. Thus the appropriate 
distribution of the t quarks in this chain should be the initial 
distribution of the X particles times the relevant branching 
ratio. This skews the expected v T distribution from solar X 
annihilation to higher energies. The modifications are dis- 
cussed in detail in (MacGibbon et al., 2001). 



4 Conclusions 

Previous work has assumed that the v T component in TeV 
- UHE hadronic decays is negligible compared with the v e 
and i>n components, at all neutrino energies. On examining 
QCD behaviour at accelerator energies and assuming similar 
behaviour continues to higher energies, we find this is not so. 
In particular, for neutrino energies in the decade below the 
energy of the decaying particle, the tau neutrino component 
is of similar magnitude to the v e and components. Below 
these energies the v T spectrum exhibits a slope of slightly 
less than E~ x / 2 compared with the E~ 3 ^ 2 slope of the v e 
and spectra. This analysis modifies the expected spectra 
in many astrophysical and cosmological high energy neutrino 
production scenarios. As a consequence, the observation of 
a significant v T to ratio at a given energy in high energy 
neutrino telescopes and detectors may be due to hadronic de- 
cay at the source and not — > v T oscillation in transit. 
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